Origin of the Universal Roughness Exponent of Brittle Fracture Surfaces: 
Correlated Percolation in the Damage Zone 
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We suggest that the observed large-scale universal roughness of brittle fracture surfaces is due 
to the fracture process being a correlated percolation process in a self-generated quadratic damage 
gradient. We use the quasi-static two-dimensional fuse model as a paradigm of a fracture model. 
We measure for this model, that exhibits a correlated percolation process, the correlation length 
exponent v m 1.35 and conjecture it to be equal to that of uncorrelated percolation, 4/3. We then 
show that the roughness exponent in the fuse model is £ = 2v/(l + 2v) = 8/11. This is in accordance 
with the numerical value ( = 0.75. As for three-dimensional brittle fractures, a mean-field theory 
gives v — 2, leading to £ = 4/5 in full accordance with the universally observed value £ = 0.80. 

PACS numbers: 83.80.Ab, 62.20.Mk, 81.40.Np 
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Fracture surfaces in brittle materials show surprising 
scaling properties Q. These were first seen in the mid- 
eighties Q . They manifest themselves through self-affinc 
long-range height correlations. That is, the conditional 
probability density p(x,y), i.e. the probability that the 
crack surface passes within dy of the height y at position 
x when it had height zero at x = 0, shows the invariance 



X c p(Xx,X <: y) = p(x,y) 



(1) 



where £ is the roughness exponent. In the early nineties 
increasing experimental evidence hinted at the roughness 
exponent not only existed, but had a universal value of 
about 0.80 j|. The experimental picture today is even 
more complex: a) a second, smaller roughness exponent, 
approximately equal to 0.5 has been observed on small 
length scales, with a clear crossover length between the 
two regimes j|; b) the growth of the roughness from 
an initial straight notch shows anisotropy |5J and a two- 
regime process in case of quasi-brittle material such as 
wood ||; c) materials like sandstone for which the frac- 
ture is strongly transgranular, show only a ( = 0.5 self- 
affine scaling |j. Simultaneously with these experiments, 
theoretical and numerical work have been produced at a 
steady rate with the aim of: (1) understanding why there 
is a self-amne scaling of the roughness, (2) why there 
should be universality of the roughness exponent and (3) 
how to unify observations and modelling . 

It is the aim of this letter to present a new possible ex- 
planation for the observed universal roughness of brittle 
fracture surfaces at larger scales. We present our ideas 
using a paradigm of fracture model: the quasi-static fuse 
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model Dynamical fuse models have been proposed 
and studied in the work of Sornette and Vaneste |T^] . The 
quasi-static fuse model consists of a lattice where each 
bond is an ohmic resistor as long as the electrical current 
it carries is below a threshold value. If the threshold is 
passed, the bond burns out irreversibly. The threshold t 
of each bond is drawn from an uncorrelated distribution 
pit). The lattice is placed between electrical bus bars and 
an increasing current is passed through it. Numerically, 
the Kirchhoff equations are solved with a voltage differ- 
ence between the bus bars set to unity. The ratio between 
current ij and threshold tj for each bond j is calculated 
and the bond having the largest value, max.j(ij/tj) is 
identified and subsequently irreversibly removed. 

In the limit of inifinite disorder — i. e. when the 
threshold distribution is on the verge of becoming non- 
normalizable, e. g. p(t) oc t~ a ~ x , where 1 < t < oo in 
the limit of a — ► 0, the fuse problem becomes equivalent 
to a bond percolation problem flTo| . At more narrow dis- 
orders, a rich phase diagram appears which is controlled 
by two parameters, the exponent a which controls the 
threshold distribution tail towards infinitely large thresh- 
old values and the exponent (3 which controls the tail of 
the threshold distribution towards zero: p(t) cx t~ 1+l3 
where < t < 1 jl3). For smaller values of either a 
or /3, the fuse model still shows behavior very similar to 
percolation: The lattice stops conducting after a finite 
percentage of bonds have burned out even when the lat- 
tice size is extrapolated to infinity. Close to breakdown, 
critical exponents may be defined precisely as in the per- 
colation problem. However, as the breakdown process in 
the fuse model is highly correlated, there is no reason to 
expect these exponents to be equal to those found in the 
percolation problem. At even smaller disorders, localiza- 
tion sets in. 

When the disorder is broad enough so that the fuse 
model behaves in a percolation-like like manner, there is 
a diverging correlation length £ oc \p — p c \~ u , where p is 
the density of broken bonds and p c is the density at which 
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FIG. 1: Survival probability as a function of density of broken 
bonds plotted against (p—pcJL 1 ^, where p c = 0.3735 and v — 
4/3 gives a good data collapse for different lattice sizes. The 
threshold distribution was p(t) oc t~ 1+ " on the unit interval 
where ft — 1/10. The number of samples for each lattice size 
varied from 2000 for L = 10 to 80 for L = 60. 



FIG. 2: 50% survival probability plotted against inverse lat- 
tice size to the power 0.75 for threshold distributions on the 
unit interval, p(t) oc t~ 1+/3 where (3 = 1/10 (x) and p = 1/3 
(+). The straight lines extrapolate to p c — 0.3735 and 
p c = 0.252 respectively. 



an infinite lattice breaks down. For classical percolation, 
v = 4/3 [Q. For the fuse model away from the infinite- 
disorder limit, v has not been measured. Three scenarios 
are possible for the value of v. (1) v depends on the 
disorder. Hence, it is not a universal quantity. (2) v is 
independent of the disorder but is different from 4/3. In 
this case, the fuse model defines a new universality class 
different from standard percolation. (3) v is the same in 
the fuse model as in standard percolation. Thus, the fuse 
model is in the universality class of percolation. 

In order to determine which of these three scenarios is 
correct for the two-dimensional fuse model, we studied 
the survival probability of lattices for different system 
sizes and different disorders. In Fig. |], we show survival 
probability for the threshold distribution p(t) oc t~ 1+l3 
when < t < 1, where f3 = 1/10 as a function of density 
of broken bonds for different lattice sizes. The collapse of 
the curves obtained for differents sizes shows both that 
the survival probability is converging on a step function 
at a finite p = p c , and that an estimate of the coefficient 
1/v = 0.75. Indeed, we expect the survival probability 
to scale as i" 1 /". In Fig. |, we confirm the estimate of 
v by showing the 50% survival probability, as a function 
of lattice size for this threshold distribution and for the 
threshold distribution p{t) oc t~ 1+l3 on the unit interval, 
where (3 = 1/3. Finite size scaling dictates that the ef- 
fective density at which 50% of the lattices survive, p s , 
behaves as 

Ps=Pc~ j^; , (2) 

where c is a constant. By adjusting the value of v until 
a straight line ensues, we determine the value of v. We 



find v = 4/3 fits the data very well. These results are 
consistent with scenario (3) above: The two-dimensional 
fuse model is in the same universality class as classical 
two-dimensional percolation. 

If one tries to determine the scaling properties of the 
final crack for disorders that are so broad that the system 
behaves as regular percolation, the roughness exponent 
will be one, since the final crack will be fractal with no 
anisotropy and so the width of the crack will essentially 
be that of the lattice itself. However, with more narrow 
threshold distributions, a non-linear gradient develops in 
the damage profile in the average current direction. That 
is, if y is the average current direction (which is normal 
to the two bus bars), then damage density averaged in 
the orthogonal x direction, (p) (y) takes the form 

( P ){y)=P S -A[y^j 2 , (3) 

where A is a positive constant that depends on the width 
of the threshold distribution and l y is the width of the 
damage distribution. The damage profile must surely be 
quadratic as the system must be statistically mirror sym- 
metric about y c , where the maximum damage occurs. At 
breakdown, the maximum damage pf is equal to the crit- 
ical damage density p c and can be expressed in terms of 
the correlation length £: ((p) — p c ) oc . As proposed 
by Sapoval et al. for percolation in a gradient pjfl , we 
suggest to consider the region along the damage zone that 
is at a distance corresponding to the correlation length: 
(y — y c ) oc £. Stating that the crack roughness is propor- 
tional to the correlation length: £ w w and solving Eq. 
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FIG. 3: Damage profile (p)(y) normalized so that its its max- 
imum is set to unity plotted against y/L for different lattice 
sizes L and for the threshold distribution p(t) oc t~ 1+l3 on the 
unit interval and where j3 = 1. The curve is a quadratic best 
fit based on the L — 32 data in accordance with Eq. (]3j). 



FIG. 4: Width of the damage distributions shown in Fig. H 
(+) and one based on the threshold distribution p(t) oc t~ 1+ ~ 
on the unit interval, whith j3 = 1/3 (o) plotted against L. 
The straight lines are linear fits to the data. 



(||) with the above conditions yield 



w ~ ^"AW 



(4) 



The width of the damage profile, l y , must be propor- 
tional to the length of the system L. The reason for 
this is that each broken bond creates a disturbance in 
the average current field that enhances the probability 
for a new bond to break in a finite-width cone which 
stretches out from each side of the bond in the direction 
approximately orthogonal to the average current direc- 
tion. Hence, as long as the current enhancement is not 
sufficient to induce crack coalescence and create an un- 
stable crack tip, the damage zone will spread in the new 
cones in a random fashion. This leads to 



l y oc L 



(5) 

In Fig. [|, we show the damage profile averaged over 
many samples and for many lattice sizes plotted against 
y/L. We note from Fig. ^that the profiles clearly follow 
Eq. (||). The collapse of the damage profiles shows that 
they are functions of the combination y/L and L docs 
not enter in any other way. This results is confirmed in 
Fig. [| where the width of the damage zone l y is plotted 
versus the system size L for two different threshold distri- 
butions. Both show a good linear behavior in accordance 
with Eq. (Hj. Hence, the width of the crack scales as 



L 



2v/{l+2v) 



(6) 



We therefore conclude that the fracture rougness expo- 
nent of the fuse model is 



c 



2v 



\ + 2v 



11 



0.73 , 



(7) 



where we have on the right hand side assumed that v — 
4/3, the standard percolation value. In Ref. g, ( was 
measured to be about 0.75 in the two-dimensional fuse 
model. Hence, there is very good agreement with Eq. (Q). 

We now extend our argument to the general case of 
quasi-static 3D brittle fractures in heterogeneous mate- 
rials. No measurements of v exist for the brittle fracture 
problem. However, a recent mean-field theory obtained: 
v = 2 |]l7| , making v very different from the value found 
in standard three-dimensional percolation, v = 0.88. Us- 
ing Eq. (0), we arrive at 



c 



(8) 



which is indeed in excellent agreement with the exper- 
imentally observed roughness exponent for large scales, 
C = 0.80. 

Why should this theory be applicable only to the large- 
scale exponent observed in brittle fracture, and not the 
exponent seen at small scales? It is the large-scale expo- 
nent that describes the correlated behavior of the damage 
field which finally will lead to the large-scale properties of 
fracture surface. The smaller exponent, which is close to 
0.5, describes opening of small cracks, and may be caused 
by corrugation waves propagating elastically along the 
crack front fpTsj . At larger scales, where a roughness ex- 
ponent equal to 0.8 is observed, these waves are too weak 
to influence the system. At this larger scale, we propose 
that it is a correlated gradient percolation process which 
is responsible for the value of the roughness exponent. 

This work was partially funded by the CNRS PICS 
contract #753 and the Norwegian Research Council, 
NFR. We thank Fernando Oliveira and the ICCMP for 
financial support during our stay in Brasilia. Discussions 



4 



with K. J. Mal0y and G. G. Batrouni are gratefully ac- knowledged. 



[1] 
[2] 



[3] 



[4] 

[5] 
[6] 
[7] 

[8] 
[9] 



E. Bouchaud, J. Phys. Condens. Matt. 9, 4319 (1997). 
B. B. Mandelbrot, D. E. Passoja, and A. J. Paullay, Na- 
ture, 308, 721 (1984); S. R. Brown and C. H. Scholz, J. 
Geophys. Res. 90, 12575 (1985). 

E. Bouchaud, G. Lapasset, and J. Planes, Europhys. 
Lett. 13, 73 (1990); K. J. Mal0y, A. Hansen, E. L. Hin- 
richsen, and S. Roux, Phys. Rev. Lett. 68, 213 (1992); 
J. Schmittbuhl, S. Gentier, and S. Roux, Geophys. Res. 
Lett. 20, 639 (1993); B. L. Cox and J. S. Y. Wang, Frac- 
tals, 1, 87 (1993); J. Schmittbuhl, F. Schmitt, and C. H. [10] 
Scholz, J. Geophys. Res. 100, 5953 (1995). 

P. Daugier, S. Henaux, E. Bouchaud, F. Creuzet, Phys. [11] 
Rev. E 53, 5637 (1996); P. Daugier, B. Nghiem, E. 
Bouchaud and F. Creuzet, Phys. Rev. Lett. 78, 1062 
(1997). [12] 
J. Schmittbuhl, S. Roux and Y. Berthaud, Europhys. 
Lett. 28, 585 (1994). 

S. Morel and J. Schmittbuhl and J. Lopez and G. [13] 
Valentin, Phys. Rev. E 58, 6999 (1998). 

J. M. Boffa, C. Allain and J. P. Hulin, European Phys. [14] 
J. - App. Phys. 2, 281 (1998); Y. Meheust, PhD Thesis, 
Universite Paris XI, (2002) [15] 
A. Hansen, E. L. Hinrichsen and S. Roux, Phys. Rev. 
Lett. 66, 2476 (1991). [16] 
J. P. Bouchaud, E. Bouchaud, G. Lapasset and J. Planes, 
Phys. Rev. Lett. 71, 2240 (1993); A. B. Mosolov, Eu- [17] 
rophys. Lett. 24, 673 (1993); E. Bouchaud and J. P. 
Bouchaud, Phys. Rev. B 50, 17752 (1994); S. Ra- 
manathan, D. Ertas and D. S. Fisher, Phys. Rev. Lett. [18] 
79, 873 (1997); G. G. Batrouni and A. Hansen, Phys. 
Rev. Lett. 80, 325 (1998); V. I. Raisanen, E. T. Seppala, 



M. J. Alava and P. M. Duxbury, Phys. Rev. Lett. 80, 329 

(1998) ; V. I. Raisanen, M. J. Alava and R. M. Nieminen, 
Phys. Rev. B 58, 14288 (1998); C. W. Lung, J. Jiang, 
E. K. Tian and C. H. Zhang, Phys. Rev. E 60, 5121 

(1999) ; E. T. Seppala, V. I. Raisanen and M. J. Alava, 
Phys. Rev. E 61, 6312 (2000); B. Skjetne, T. Helle and A. 
Hansen, Phys. Rev. Lett. 87, 125503 (2001); F. Barra, H. 

G. Hentchel, A. Levermann and I. Procaccia, Phys. Rev. 
E, 65, 045101 (2002). 

S. Roux, A. Hansen, H J. Herrmann and E. Guyon, J. 
Stat. Phys. 52, 237 (1988). 

H. J. Herrmann and S. Roux, Statistical Models for 
the Fracture of Disordered Media (Elsevier, Amsterdam, 
1990). 

D. Sornette and C. Vanneste, Phys. Rev. Lett. 68, 612 
(1992); L. Lamaignere, F. Carmona and D. Sornette, 
Physica A 241, 328 (1997). 

A. Hansen, E. L. Hinrichsen and S. Roux, Phys. Rev. B 
43, 665 (1991). 

D. Stauffer and A. Aharony, Introduction to percolation 
theory (Francis and Taylor, London, 1992). 

B. Sapoval, M. Rosso and J.F. Gouyet, J. Phys. Lett. 
(France) 46, L149 (1985). 

S. Roux, A. Hansen and E. L. Hinrichsen, J. Phys. A 23, 
L1253 (1990). 

R. Toussaint and S. Pride, Fracture of disorderd solids 
in compression as a critical phenomenon: III Analysis of 
the localization transition, Phys. Rev. E (submitted). 

E. Bouchaud, J. P. Bou chaud. D. S. Fish er. S. Ra- 



manathan and J. R. Rice, |sond-mat/010826l| (2002) 



